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Abstract 

The task under investigation is determination of disposition 
of emission source according to the temporal delays of the 
signals accepted in the network of sensors diversity. The 
criterion of efficiency for optimization of placing of sensors 
and the standard examples of network configuration are 
worked out. 
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Introduction 

Source localization by means of passive location 
systems is one of the central problems in sonar, 
navigation, wireless communication, geophysics and 
many others. Thus, important attention is paid to a 
choice of criteria and methods of combining the data 
for increasing the radiating sources location efficiency, 
taking into account a spatial configuration of a sensors 
network [1-3]. 

Nowadays the algorithms of radiation sources 
coordinates determination based on measurement of 
temporal time delays between signals, arriving on 
spatially delivered sensors are widely applied [1-5]. In 
case of remote sources of radiation it is supposed that 
the wave front is close to plane and coordinates are 
calculated by triangulation methods [2, 3]. However in 
conditions when the extent of basis is commensurable 
with the distance to a radiation source, the wave front 
is spherical. It leads to creation of the mathematical 
model, that is containing system of equations, 
nonlinearity of which generates several decisions and 
ambiguity in case of determination of the radiation 
source position [4-5]. 


Simulation of the Sensor Network 

The present investigation is devoted to research of 
influence of a spatial configuration of a sensor 
network on the accuracy of signal radiation source 
coordinates calculation in the controlled zone. 

It is supposed, that false peak signals are absent, 
mutually correlative functions of processes on outputs 
of the selected couples of sensors calculated correctly, 

with them the temporal signal delays T..,T km 
between pairs of sensors are calculated i and j , k 
and Jfl . Calculation of coordinates (x,jp) of a source 
of radiation which is carried out on the basis of the 
decision of system of the non-linear (hyperbolic) 
equations which have been worked out on temporal 

time delays of found signals from L pairs of sensors 
is necessary. 

The most known model of system of direction finding 
for determination of coordinates of the source of 
radiation consists of the L =2 bases containing 2 
sensors everyone [2]. The decision (x, y) in Cartesian 
coordinates about location of the source of radiation in 
case of arbitrary known coordinates of sensors 

X i ,y i )i=l,.-A is defined from relation of temporal 
time delays T l2 ,T 34 of each of couples, and distances 
r x , r 2 , r 3 , r 4 from a radiation point to 
sensors r 34 = r 3 - r 4 = c T 34 , r l2 = r x - r 2 = cT l 2 
where C - the speed of a sound wave. Coordinates 
(x, y) of object of radiation are described by a 
mathematical model in the form of system of two 
hyperbolic equations [6, 7] which have been worked 
out according to the Pythagorean theorem 
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r n = r r r 2 = ^(x, -x) 2 +(j, -j) 2 - V(* 2 -*) 2 + 0> 2 - j) 2 

/ 

r 34 =r 3 -r 4 = j(x 3 -x ) 2 + (y 3 -y ) 2 = ^J(x 4 -x ) 2 + (y 4 - y ) 2 . ^ 

For the method of computation of coordinates on time- 
differences-of-arrival (TDOA) some types of the 
ambiguities exist arising from the following reasons. 

• Difference of the discrete space in polar and 
Cartesian coordinate systems. Thus, one couple 
temporal time delays t = T 12 ,t = T 34 in (1) 

corresponds to several nearby to the located 
elements of permission (cells) 

t,T — > X , V / = 1,..,7V in Cartesian coordinate 

7 i 7 y i 7 77 t,r 

system. 

• Direct decision of system (1) method of squaring. 
Thus there are some decisions generating speckled 
false points, significantly remote from the 
coordinate of a source of radiation. In particular, 
essential ambiguity takes place in case of the 
analytical decision (1) as a result of double 
exponentation of the equations in a square. It 
generates to 4 decisions provided by vectors 
S. = (x i ,y i ) T , i = l v .,4, from which 3 false and 
don’t correspond to the true coordinates of a source 
of radiation. It is necessary to mark that receiving 
exact analytical expressions even for the simplest 
implementation of a network from L =2 bases 
containing randomly located M =3 sensors (two 
couples sensors, one of them the general), results in 
ambiguity in the form of two decisions [8]. 

For simplification of calculations of coordinates with a 
margin error less than a half of the size of an element 
of permission approach which is based on grid 
sampling of a card of compliance of temporal time 
delays T x 2 ,T 3 4 — > x , y to each point in the Cartesian 

coordinate system calculated in advance is offered (1). 
It allows to exclude the need of the decision of system 
of the hyperbolic equations if the size of a grid of 
coordinates of the purposes for all possible time delays 

T l2 J 34 taking into account the size of a cell of a 

controlled zone is limited. It is easy to define the 
requirements to the size of a grid and to the memory 
necessary for its storage. An example of a big zone of 
monitoring, characteristic for the acoustic 
reconnaissance [2, 9] which is carried out on a square 
D = D "10 km is reviewed. In case of sampling rate of 

f = 1000 Hz and the speed of a sound wave of v=330 
m/s, the size of an element of permission on range 
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make d — vlf = 0,33 m. Then storage of cells 
N N =D D / d 2 & 10 9 °f a card of compliance (the 

cell occupies 8*4=36 bytes) requires about 36 Gbytes of 
memory that is implemented on the modern element 
basis. 

For research of ambiguity an example of nonoptimal 
placement of the network of sensors in a controlled 
zone (Figure 1) in the form of a square with the side is 

reviewed X = V =2 m. Bases are located 

max y max 

parallely, and steam 1 and 2 of sensors make the first 
basis of the size 5, = 5(1,2) = -x 2 ) 2 +(y t -j 2 ) 2 , 

and steam 3 and 4 of sensors - the second basis 
B 2 = 5(3,4) = ^(x^x^ 1 +(y 3 -y 4 f • The audible 
signal is digitized with sampling rate of f =1000 Hz 

that at a speed of distribution of a sound wave of 330 
m/s corresponds to a permission element on range of 
d = v/ / =0,33 m. 

From Figure la it is easy to note that any line item of a 
source of radiation on the shaped line passing through 
centers of bases, leads to identical time delays of 
signals T l2 -t x ~t 2 - 0, T 3 4 =t 3 ~t 4 = 0 (the number 
of elements is determined by a scale of intensity). Then 
the set of decisions (1) represents set of all elements 
located on the shaped line. It generates essential 
ambiguity in position fix of a source of radiation as the 
source is only in one cell on this straight line, and 
remaining decisions are false. 

In Figure lb the card of a zone of monitoring with 
location of ambiguous elements in Cartesian 
coordinate system on which sampling of elements of 
permission is carried out is provided. Multiplicity of 
ambiguity is determined by a scale of intensity 
(number of the false decisions corresponding to one 
couple time delays). Apparently, in case of a time 

delay t = ZJ 2 =0 on basis B { and T = T 3 4 =0 on basis 
B 2 the greatest multiplicity of ambiguity, to equal 
N t = N Q 0 = 13 elements of permission or the size of 

a controlled zone takes place. In case of other time 
delays ambiguity is less, and on white squares of 
coordinate of a source are defined by the single 
decision. From the analysis of Figure la and Figure lb 
follows that it is necessary to define the configuration 
of sensors providing the minimum number of 
elements of permission, time delays corresponding to 
everyone couple. 

For comparing of configurations of systems it is 
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necessary to define the criterion of ambiguity 
considering precision characteristics of an assessment 
of location of the purpose. 




FIG. 1 EXAMPLE OF AMBIGUOUS DETERMINATION OF 
COORDINATES: CONFIGURATION OF SENSORS AND A 

RADIATION SOURCE ON THE PLANE OF A ZONE OF 
MONITORING (A); LOCATION OF AMBIGUOUSLY DEFINED 
ELEMENTS ON THE PLANE OF A ZONE OFMONITORING (B) 

If the decision on signal existence on temporal time 
delays t = T 12 ,t = T 34 is made, in Cartesian 

coordinate system it can correspond N t T > 0 to 

several points S t (x i ,y i ) 9 i = • Mean value 

and dispersion of an assessment of coordinates are 
defined by expressions: 






( 2 ) 


In case of the single decision parameters of delay 
correspond to t,T —>x,y one element of permission 

(cell) the minimum error CT t T =0 is reached. With 
ambiguity when to one couple temporal time delays 

t,T —> X . , y . , i = 1 9 .. 9 N tT corresponds two and more 
close cells in Cartesian coordinate system, takes place 


C t >0, and value of ambiguity is equivalent to 

lowering of resolution capability and accuracy of an 
assessment of location of the purpose. If cells don’t 
adjoin, mean squared error (MSE) CF ^ significantly 

increases that corresponds to appearance of false 
targets. 


These features allow to evaluate numerically 
ambiguity and to define efficiency of placement of 
sensors, as the amount of the errors enveloping 
possible provisions of the purpose. As criterion of 
ambiguity the error square, equal to the amount of 
dispersions on all possible time delays in a controlled 
zone is accepted: 



Z L max 
T—l 




(3) 


where W tz - a weight multiplier of importance of the 
cells enveloped by time delays T ; 
T — B / cA — L2 - the maximum time delays of 
signals on bases B. i - pairs of sensors. 

Then determination of a configuration M of sensors 
is brought together to the choice of coordinates of 
sensors in the form of vectors X \Y of the size M 

and to the choice of bases B n i - 1 ,L on pairs of 

sensors. The task of optimization is brought together 
to functionality minimization from MSE of an error 
and looks like 

miner (X,Y,B l9 B 2 ,..,B L ) (4) 

in case of restrictions B < B i = 1 L' 

i max ’ 5’ 

x <x , v <v ,i = LM on the maximum 

i max 5 si J max 5 ? 

sizes: bases B and controlled zone r v 

i max 5 J max 

The optimum is reached as a result of integer search of 
all values X,Y and bases B . , i = \ L on a grid of a 
controlled zone with a permission element d =v/ f • 

Generally, the maximum possible quantity of bases is 
defined by number of sensors L=M(M—Y)/2, 

and the mathematical model of system of direction 
finding (1) becomes complicated and contains to the 
L equations 

r tJ = r- r , = cT tj =J(x, -x ) 2 +(y, - j ) 2 -^(x, -x ) 2 +(j. - j ) 2 

/ 

i,j = l,2,-,M; i*j- (5) 
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In practice the situation in which the configuration is 
set often takes place, and it is necessary to select pairs 
of sensors for bases in case of which ambiguity 
minimum. 

Let’s begin efficiency research on the configuration 
containing the minimum quantity of sensors M= 3, 
necessary for purpose position fix. In Figure 2a the 
configuration "angle" on the sides of a square zone in 
which two bases are applied is shown, first of which 
contains 1 and the 2nd sensors B x = 5(1,2), the second 
-2 and the 3rd sensors B 2 = 5(2,3) • 

Apparently, zones of ambiguity contain to 3 elements 
of permission and concentrate in 3 zones. Increase of 
efficiency is possible at the expense of use of all 
possible bases on this configuration of sensors. The 
result of application of the third additional basis for 
which it is necessary to enter the third equation in (1) 
is given in Figure 2b. It significantly reduces 
ambiguity and the MSE decreases by 4 times. In Figure 
2c ambiguity zones for system with the high 
resolution, different system in Figure 4 sampling rate 
increase by 5 times are provided. From their 
comparing follows that the same circuits of ambiguity 
that allows to conduct researches on models with low 
resolution and, thereby, to reduce computing costs of 
search of an optimum configuration remain. 

In Figure 3 examples of placement of 4 sensors on a 
configuration "cross" are given 5 — 5(1,3) , 

B 2 = 5(2,4)- For the system provided in Figure 3a, 

bases are square diagonals. In this case the extent of 
basis exceeds the side of a square of the zone of 
monitoring that allows to exclude completely 
ambiguity and C =0. If sensors are distributed on the 
circle inscribed in a square (Figure 3b), the extent of 
basis decreases to the side of a square and ambiguity 
at the edges of a monitoring zone takes place. The 
essential increase in ambiguity arises in case of small 
bases of rather controlled zone (Figure 3c). In practice 
such situation corresponds to use of systems of 
detection of snipers for which microphones are placed 
on a helmet of the security guard with bases about 0,2 
meters [10], and the size of a zone of monitoring 
occupies hundreds meters. From Figure 3c follows that 
the worst case situation corresponds to location of the 
shooter on square diagonals, only the direction of a 



FIG. 2 EXAMPLES OF PLACEMENTS OF SENSORS AND 
AMBIGUITY ZONES FOR CONFIGURATIONS OF THE BASES 
MADE ON 3 SENSORS: APPLICATION OF TWO BASES (A); 
THREE BASES (B) ; THREE BASES AND THE INCREASED 
SAMPLING RATE OF A SIGNAL (C) 


ambiquity is absent 

a =0 

optimum placement 
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FIG. 3. EXAMPLES OF PLACEMENTS OF SENSORS AND 
AMBIGUITY ZONES FOR CONFIGURATIONS OF THE BASES 
MADE ON 4 SENSORS: A CONFIGURATION "CROSS" 
APPLICATION OF TWO BASES (A);SENSORS ARE 
DETRIBUTED ON A CIRCLE (B); CONFIGURATION WITH THE 
REDUCED BASES (C) 

Example of other placement of 4 sensors is the 
configuration "line" with the even stride which 
optimization is illustrated in Figure 4. Apparently, the 
case of use of L =2 two bases = 5( 1,2) and 

B 2 = 5( 3,4)/ provided in Figure 4a, possesses several 

shortcomings connected to the small extent with the 
bases and their nonoptimal placement in the zone of 
monitoring. 




FIG. 4 EXAMPLES OF CONFIGURATIONS OF THE BASES 
LOCATED ON THE LINE: TWO BASES (A); THREE BASES (B); 
FOUR BASES (C) 

Optimization at the expense of a choice of bases 
without change of position of sensors is performed. 

Application of the third basis (Figure 4b) L =2 look 
allows B 3 = 5( 2,3) f° reduce an error by 4 times, the 

fourth basis (Figure 4c) L= 4 5 4 = 5(1,4) - by 1,5 
times. It is necessary to mark that for big zones of 
monitoring, for example, in tasks of artillery acoustic 
reconnaissance [2] with zones to 10-30 km, application 
of big bases of a look isn't always possible 
B 4 = 5(1,4). It is connected with possible lowering of 
correlation of signals because of distortions of the 
signals entered by different propagation paths to 
sensors, and also to relation reduction signal/noise on 
inputs of sensors. 


Conclusions 

The conducted researches on calculations of location 
of radiation source allow to draw the following 
outputs: 

• The criterion is developed and the task of 
optimization of the sensor network configuration 
taking into account restrictions on the sizes of a 
controlled zone is formulated; 
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• It is set that ambiguity minimum in case of 
execution of the following requirements to 
arrangement of sensors: application of the greatest 
possible number of bases and their sizes, support 
of orthogonality of lines of bases by means of 
approximation to a configuration "cross". 

The considered method of mathematical simulation of 
a configuration allows to realize reasonably creation of 
a sensor network and to predict precision 
characteristics of a passive location in a controlled zone. 
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